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ABSTRACT 

Hierarchical triple-star systems are expected to form frequently via close binary-binary encounters 
in the dense cores of globular clusters. In a sufficiently inclined triple, gravitational interactions be¬ 
tween the inner and outer binary can cause large-amplitude oscillations in the eccentricity of the 
inner orbit (“Lidov-Kozai cycles”), which can lead to a collision and merger of the two inner compo¬ 
nents. In this paper we use Monte Carlo models of dense star clusters to identify all triple systems 
formed dynamically and we compute their evolution using a highly accurate three-body integrator 
which incorporates relativistic and tidal effects. We find that a large fraction of these triples evolve 
through a non-secular dynamical phase which can drive the inner binary to higher eccentricities than 
predicted by the standard secular perturbation theory (even including octupole-order terms). We 
place constraints on the importance of Lidov-Kozai-induced mergers for producing: (i) gravitational 
wave sources detectable by Advanced LIGO (aLIGO), for triples with an inner pair of stellar black 
holes; and (ii) blue straggler stars, for triples with main-sequence-star components. We find a realistic 
aLIGO detection rate of black hole mergers due to the Lidov-Kozai mechanism of ^ 1 yr“^, with 
about 20% of these having a hnite eccentricity when they first chirp into the aLIGO frequency band. 
While rare, these events are likely to dominate among eccentric compact object inspirals that are 
potentially detectable by aLIGO. For blue stragglers, we find that the Lidov-Kozai mechanism can 
contribute up to ^ 10% of their total numbers in globular clusters. 


1. INTRODUCTION 

Observations indicate that a large fracti on of stars 
reside in triple systems (jRiddle et al.ll201^ and refer¬ 
ences therein). In a dynamically stable non-coplanar 
triple the gravitational interaction between the inner 
and outer binary drives periodic/quasi-periodic varia¬ 
tions of the mutual inclination between th e two orbits 
and of the e ccentricity of the inner binary (jLidovlFTOBlI : 
iKozail 11962). Such ecc entricity oscillations, first dis- 
covered bv iLidovI (I1961D . can lead to a close approach 
between the inner binary components which can pro¬ 
duce a variety of exotic astrophysical phenomena, in- 
cluding compact-object and stellar binary mergers (e.g.. 

Ford et al.l2004lMiller fc Hai:niltonll2002t I Antonini et al'l 

20irt iNaoz fc Fabrvckvl I2014D . X-rav binaries, type la 
supernovae and G amma-ray bursts (iThompsonl 120111 : 
iHamers et al.l[2M^ . 

Lidov-Kozai (LK) oscillations are usually described as 
a secular phenomenon, occurring on timescales much 
longer th an the orbital period of the outer binary or¬ 
bit (e.g.. iFord et al.l 12004 iFabrvckv fc Tremaine! I2007t 
INaoz et al.ll201.‘I[ ). The standard approach is to make use 
of the so-called orbit-averaging technique in which one 
computes the mutual gravitational interaction (torques) 
between two mass weighted ellipses instead of point 
masses on orbits. The basic assumption is that the evolu¬ 
tion of the orbital parameters occur on a timescale much 
longer than the dynamical timescales of the system (e.g., 

The orbit averaged equations of motion can describe 
the evolution of a wide range of systems, but be- 


come inaccurate once th e triple hierarchy is moderate. 
lAntonini fc PeretsI (|2012D showed that the orbit-average 
treatment, even at the octupole level of approximation, 
becomes inaccurate in hierarchical systems where the ter¬ 
tiary comes closer to the inner binary than a certain 
distance. At the high eccentricity phase of a LK cycle 
the timescale for the angular momentum perturbation 
due to the tertiary can become comparable to or shorter 
than the dynamical timescales of th e system rendering 
the s e cular approach no longer va l id dAntonini fc Perets 
2012: iKatz fc P ond 20l3 ISetol [201,4 lAntoiiini et al. 


20I4 : IAntognini et al. 12014 iBode fc Wesfell20I4[) . In this 


case, the binary can be driven to a much closer sepa¬ 
ration than the secular theory would otherwise predict. 
This can have important consequenc es for the evolution 
of the inner binary . F or example. lA ntonini fc PeretsI 
(j2012fl . iSetol (j201.‘lf ) and lAntonini et al.l (j20l4 ) demon- 
strated that in this non-secular regime, a compact ob¬ 
ject binary can be pushed inside the 10 Hz gravitational 
wave (GW) frequency band of Advanced Ligo (aLIGO) 
detectors with a finite eccentricity. Similarly, in the case 
of stellar bin aries, the two stars could be driven to a di - 
rect collision. INaoz fc Fabrvckvl (j20l4 ) and lPereii^ (j2015[ ) 
suggested that direct stellar collisions driven by the LK 
mechanism might have an important role in the forma¬ 
tion of blue straggler stars (BSSs) in stellar clusters. 

Globular clusters (GCs) are a natural environment 
for forming triples and higher order multiples. On one 
hand, the dense stellar environment of GCs is hostile to 
triples. In fact, given the high densities of GCs, triples 
may be involved in close encounters which tend to de- 
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stroy them. On the other hand, the high stellar densi¬ 
ties that characterize the core of GCs favor dynamical 
binary-binary interactio ns that naturally lead to the for - 
mation of stable triples (jSiefurdsson fc HernauistllTQQ.IIl . 
including black hole (BH) triples and mixed, compact 
object-star, triples that could not otherwise form. The 
inner binaries in the triples may be driven to merge be¬ 
fore their next interaction. Interestingly, such mergers 
have been suggested as a channel for the format i on of 
compact-object inergers (|Miller fc Hamilto'i] l2002t iWenI 
I2003I lArca-Sedd^ I2016D , and more recently as a chan¬ 
nel for the formation of BSSs through mass transfer 


and merger ( 

Perets & Fabrvckvl2009l:ILeie:h & Sillsll2011l: 

iGeller et al.l 

20131 Naoz & Fabrvckvll2014: IPeretsll2015l). 


In this paper we study the evolution of dynamically 
assembled triples inside GCs with particular focus on 
triples containing an inner main-sequence (MS) stellar 
binary and triples with an inner BH binary. We do not 
consider here the evolution of mixed star-compact ob¬ 
ject systems, or binaries containing evolved (e.g., red gi¬ 
ant) stars; we reserve a study of these latter systems to 
a future investigation. In this paper we use a Monte 
Carlo code to simulat e the long-term evolution of GCs 
(|Morscher et al.l [20151 1. and identify all triple systems 
formed during the cluster evolution. We then evolved 
these triples using a high-accuracy direct three-body in¬ 
tegrator which includes tidal and general relativity (GR) 
effects for both the inner and outer orbit. Scattering 
with other stars is also accounted for, but only in the 
sense that the triples are evolved for a time equal to the 
timescale over which they will be disrupted through grav¬ 
itational encounters with other stars. The results of the 
direct integrations are used to determine the rate of BH 
and stellar binary mergers induced by the LK mechanism 
in GCs. 

Although previous studies have shown that triples can 
form quite frequently in GCs, their detailed long-term 
dynamical evolution in the cluster environment has not 
been yet investigated. Our study sheds light on exactly 
this point, i.e. it allows us to assess for the first time the 
role of the host cluster evolution and its properties on the 
formation and dynamics of the triples. This allows us to 
place reliable constraints on the contribution of mergers 
from triples to BSS populations and on the event rate 
and properties of inspiraling BH binari es as a source of 
GW radiation for aL IGO detectors fe.g. lAasi et al.l[ 2 ?)T^ 
lAcernese et al.lU015[ j. 

We note that the dynamics of triples with either an 
inner stellar or compact object component has been in¬ 
vestigated by many authors and in a variety of contexts 
(including stellar triples, planets around stars, super- 
massive BH mergers), using the standard orbit averaged 
formalism and therefore fully neglecting non-secular ef¬ 
fects (e.g..lBlaes et al.ll2002HWenl2003HNaoz et al.ll2nilL 
I2013I iHamers et al.ll2013HNaoz fc Fabrvckvll2014H . The 
new ingredients of our study compared to previous work 
dealing with similar problems are: (i) the use of realis¬ 
tic initial conditions for our triple populations, that were 
directly obtained using Monte Carlo models of the long¬ 
term dynamical evolution of GCs; (ii) the use of accurate 
three-body integrations that allows an essentially exact 
treatment of the triple dynamics, avoiding the approxi¬ 
mations that are made when using the secular equations 
of motion. 


The paper is organized as follows. We begin with a dis¬ 
cussion of dynamically assembled triples in GCs and their 
properties in Section[2l and show that for most triples the 
orbit average approximation is expected to break down. 
In Section|3]we discuss our GC models, triple initial con¬ 
ditions and adopted numerical methods. We present our 
results on BH-triples and associated aLIGO event rates 
in Section ID In Section [5] we describe the results con¬ 
cerning stellar triples and their implications for BSS for¬ 
mation in GCs. We further discuss the implications of 
our results in Section |6l Section [3 sums up. 


2. BREAK-DOWN OF THE SECULAR APPROXIMATION 

We consider a binary of components of masses niQ and 
mi orbiting a third body of mass m 2 . We indicate the 
eccentricities of the inner and outer orbits, respectively, 
as Cl and 62 , and semi-major axes oi and 02 ; we define wi 
as the argument of periapsis of the inner binary relative 
to the line of the descending node, UJ 2 as the argument 
of periapsis of the outer orbit, and I as the orbital incli¬ 
nation of the inner orbit relative to the outer orbit. 

In what follows we only conside r triple systems that 
satisfy the stability criterion of iMardling fc AarsethI 

(Mf) . 


02 3.3 ^ f ^ \ \ 1 + 62 

01 ^ 1-62 Mf,y (1 _ g2)l/2 

x(l-0.3//7r) , 


( 1 ) 


with Mb = mo -I- mi. 

If the system satisfies the stability criterion © it is 
also reasonable to describe the dynamics of the entire sys¬ 
tem as the interaction between an inner binary of point 
masses mo and mi and an external binary of masses Mb 
and m 2 . We define the angular momenta Li and L 2 of 
the inner and outer binary: 


Li = momi 


Gai (1 



Mb 


( 2 ) 


and 


L 2 


Mb m 2 


Ga2 (1 — el) 

Mb -I- m2 


( 3 ) 


with G the gravitational constant. The dimensionless an¬ 
gular momenta are ii = Li/Li^c and ^2 = L 2 /T 2 ,c with 
Li,c = Li!— el the angular momentum of a circular 
orbit with the same semi-major axis, a^. 

If the changes in the orbital properties of the inner 
binary caused by its gravitational interaction with the 
third body occur on a timescale longer than both the 
inner binary orbital period. Pi, and the tertiary or¬ 
bital period, P 2 , it is convenient to average the equa¬ 
tions of motion over the rapidly varying mean anoma¬ 
lies o f the inner and outer orbits (e.g., iHamers et al.l 
12015(1 . The resulting double averaged Hamiltonian is 
% = kW with k = ?)Gmom^m 2 all\iMba%{\ — 
and (jLidov &: Ziglinl[l976HMiller fc: Hamiltoijr2002li 

W(wi,ei) = -2(1 - e\) -|- (1 - ei)cos^/ (4) 

-|-5eiSin^a;i(cos^/ — 1). 
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The quadrupole-level secular perturbation equations 
can be easily derived from the conserved Hamiltonian (|^. 
The resulting evolution equation of the inner binary or¬ 
bital angular momentum due to th e torque from the ter¬ 
tiary is (e.g., iLidov fc ZiglinlUOT^ : 


dt 


IStt 1712 a? ef 

“A4ai(l-ei)3/2 A 


sin^ / sin 2wi . 


(5) 


The cha racteristic timescale for the eccentricity oscilla¬ 
tions is (iHolman et al.l 119971) : 


We call this regime “moderately” non-secular. In this 
regime systems can experience eccentricity oscillations 
over the outer binary period and can be driven to 
higher eccentricities than secular theory would oth¬ 
erwise predict, which in the case of compact object 
binari es could result in su bstantially reduced merger 
times (|Antognini et al.ll2014ll . 

If gin = Siin/idiss ^ 1, the inner binary will undergo 
eccentricity oscillations on the timescale of the inner or¬ 
bital period; this can be achieved if the tertiary comes 
closer to the inner binary than the separation 


\m2j \aij 


02 ^ 2.5 

Oi I — 62 




( 10 ) 


Using Equation m and taking the relevant limit ei —^ 
I, we find that the maximal change in the inner binary 
angular momentum due to the interaction with the third 
object and over the inner binary orbital period is 


Siin 


_ 

M. ai (1 - • 


(7) 


and the change in the binary angular momentum over 
one orbit of the tertiary is 

Siout ~ (8) 

thus, the change experienced over the timescale associ¬ 
ated with the period of the outer orbit is naturally P 2 /Pi 
larger than the change experienced over the inner binary 
orbital period. 

Equation (jd]) describes the gravitational interaction of 
two weighted ellipses rather than point masses on orbits 
and is accurate as long as Wout/f'i ^ 1, be., changes of 
the inner binary angular momentum occur on a timescale 
longer than the tertiary orbital period. If instead the 
change in the inner binary angular momentum over the 
inner and outer orbital periods is of order £ 1 , then the 
system no longer meets the conditions required for the 
orbit averaged approximation, and more accurate direct 
three-body integrations are needed. 

Eollowing the discussion above, in order to address the 
reliability of the orbit averaged approach in describing 
the evolution of a triple system it is useful to compare 
the change in angular momentum of the binary to the 
critical angular momentum £diss ~ i/^Ddiss/ai^ where 
Ddiss is the scale at which other dynamical processes will 
dominate the evolution of the inner binary. Since in this 
paper we focus on triples in which the inner binary can 
be driven to a strong close interaction or even a merger, 
we identify these processes with energy loss due to GW 
radiation and tidal dissipation in the case of BH and 
stellar binaries, respectively. 

The condition that the inner binary undergoes rapid 
eccentricity oscillations over the outer orbital period can 
be written as gout = ^£out/£diss ^ 1; this requirement 
translates into a condition on the outer perturber sepa¬ 
ration. 


We call the regime defined by this equation “strongly” 
non-secul y (see alsolAntqnini fc PeretsI (|20I2D and equa¬ 
tion 7 in iKatz fc Pond (12012m . In this situation the 
outer perturber is close enough to significantly change 
the angular momentum of the eccentric binary at the last 
apoapsis passage leading to a jump in angular momen¬ 
tum of order £ 1 . The angular momentum of an eccen¬ 
tric inner binary can jump from several times £diss (far 
enough to avoid GR or tidal effects) to an arbitrarily 
small value. Equation m defines the region of param¬ 
eter space within which a compact object binary might 
be expected to ent er the aLIGO frequency band with 
a finite eccentricity (lAntonini fc Peretdl20ia ISetoll201^ 

lAntonini et al.l[2014f) or, in the case of inner stellar com¬ 

ponents, the binary members might expe rience a close 
interaction leading to a stellar collision (|Katz fc Pond 
[2nTaiProdan et al.lIM^ . 

Here we quantify the role of non-secular effects for 
triple systems forming in stellar clusters by analyzing the 
properties of all dynamically stable triple s that were pro¬ 
duced in the GG Monte Garlo models of iMorscher et all 
(I 2 OI 5 II (see Section The Monte Garlo models give 
the masses, stellar radii, semi-major axes, eccentricities 
and stellar type of all triples formed during the cluster 
dynamical evolution; these properties were used to make 
Figured] which shows the period distributions as well as 
the distribution of the ratio 02(1 —e 2 )/ai of all BH triples 
(left panels) and stellar triples (right panels) that formed 
in these models Q. 

in Figure dl we have distinguished triples which satis¬ 
fied the criterion (|9|) (green histograms) and those meet¬ 
ing the more stringent condition equation (HOI) (blue his¬ 
tograms) from the remaining systems which satisfied nei¬ 
ther conditions. In the case of stellar triples we evalu¬ 
ated equation m and (da setting Ddiss = 2(i?o + Ri) 
with Rq and Rj the physical ra dii of the inner binary 
components (IKatz &: Pond [^012f) . while for BH triples 
we ado pted a conservative d issipation scale of Udiss = 
10® cm (|Antonini et al.l[20l4) . 

Figure dl shows that for the majority of BH triples 
formed non-secular dynamical effects are expected to be¬ 
come important before GR terms can significantly affect 
the evolution of the inner binary. We conclude that the 
standard LK secular theory will fail in describing the 


^ / Mb \ / ai \ 

ai'^l-el\Mb) V Mb -I- 7712 / V -^’diss ) 

( 9 ) 


^ Hereafter we define as a stellar triple any stable triple in which 
both components of the inner binary are MS stars; BH triples are 
identified with triples in which the two inner components are both 
stellar BHs. 
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Fig. 1.— Upp er and middle panels sh ow the distribution of inner and outer orbital periods of the BH and stellar triples formed in the Monte 
Carlo models of IMorscher et ^ II2015I V The bottom panel gives the distribution of the ratio between the outer binary periapsis and inner 
binary semi-major axis. Blue and green histograms indicate systems for which the assumptions on which the orbit averaged approximation 
is based on break down before GR or tides can affect their evolution. Green histograms are systems for which gout = 5£out/^diss 7 1; blue 
histograms are systems tor which gjn = (S^in/^diss 7 1- ma-ny triples formed through dynamical interactions in star clusters, the secular 
approximation breaks down well above their typical dissipation scale, the LK equations of motion (even at the octupole order level) cannot 
accurately trace the evolution of these systems. 

of GW radiation for aLIGO detectors, or, in the case of 
stellar binaries, how the evolution might lead to a mass 
transfer event or even a stellar merger. We begin in the 
next section by describing the initial conditions and nu¬ 
merical methods adopted in our study. 

3. METHODS AND INITIAL GONDITIONS 
3.1. Numerical Methods and example cases 

The direct three-body integrations presented be¬ 
low were performed using the ARGHAIN integra¬ 
tor (|Mikkola fc Merrittl 120(1^ . ARGHAIN employs an 
algorithmically regularized chain structure and the time- 
transformed leapfrog scheme to accurately integrate the 


evolution of these systems due to the breakdown of the 
double averaging approximation. 

Only a few stellar triples in our sample satisfy the con¬ 
dition ([Toll , mainly because of their large value of Udiss- 
However a large fraction of stellar triples still met the 
condition equation ([51), meaning that the averaging pro¬ 
cedure over the outer orbit is not justified for many of 
these systems. 

In what follows we study the dynamical evolution of 
triples that form through dynamical interactions in the 
dense stellar environment of GCs. We use a high accu¬ 
racy three-body integrator to study how the inner binary 
in the triples might evolve to become a possible source 
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Fig. 2.— Dynamical evolution of a BH triple (upper panel) and 
a stellar triple (lower panel). Black lines show the results of three- 
body integrations, blue lines show the results of the octupole order 
orbit averaged equations of motion. Solid lines give the value of 
(1 — ei), dashed lines give the inner binary semi-major axis. In the 
upper panel the initial conditions for the BH triple system were: 
mo = 14.9 M 0 , mi = 13.8 Mq, m 2 = 14.2 Mq ai = 0.45 AU, 
02 = 3.9 AU, ei = 0.85, 62 = 0.51, I = 92.5°, cji = cj 2 = 7r/2. 
In the lower panel the initial conditions for the stellar triple were: 
mo = 0.7 Mq, mi = 0.5 Mq, m 2 = 2 Mq oi = 1 AU, 02 = 25 AU, 
I = 101.8°, ei = 0.1, 62 = 0.7, oji = 214.6°, a ;2 = 0. In the up¬ 
per panel the asterisks denote the minimum separation attained by 
the two BHs every orbital period. At (1 — ei) < 10“^ this system 
satisfies the condition q[^ > 1, so that the jump in periapsis over 
one inner orbit is of order the periapsis distance itself. In the lower 
panel at (1 — ei) < 0.04 the triple satisfies the condition ^out ^ 1, 
so that the averaging procedure over the outer orbit is no longer 
justified. In both examples the inner binary reaches higher eccen¬ 
tricities in the three-body simulation than in the secular integra¬ 
tion. This allows the two BHs in the upper panel to merge through 
GW energy loss, and the two stars in the lower panel to approach 
each other close enough that mass transfer would commence. The 
dot-dashed lines indicate the critical value of (1 — ei) below which 
GW radiation will dominate (upper panel), and below which the 
two stars will collide (lower panel). These two examples clearly il¬ 
lustrate the limitation of the standard secular approach even when 
including octupole order terms to the equations of motion. 


terms that account for dissipative tides as well as ap- 
sidal precession induced by tidal bulges. In order to 
do so we modified ARCHAIN including terms to the 
equations of motion which represent tidal friction and 
quadrupolar distortion. Velocity-dependent forces were 
in cluded via the generalized mid-point method described 
in iMikkola &: MerrittI (|20Q6[) : this method allows to time- 
symmetrize the algorithmic regularization leapfrog even 
when the forces depend on velocities, permitting the ef¬ 
ficient use of extrapolation methods. The tidal pertur¬ 
bation force has the form 


F = -G 


momi 


3^ 

mi 



r 

l-h3-r 

r 


( 11 ) 


where r (here set to a fixed value of 1 sec) is the time 
lag, and k (set to 0.03) is the apsidal motion constant. 

Equation (HU differs from the classical treatment of 
near-equilibrium tides to the extent that we have sup¬ 
pressed terms that are due to the difference in rotation 
rate between the binary and each component. The or¬ 
bital elements evolution equations that correspond to 
Equ a tion (fTTl) are given by Equations (9) and (10) of 
I Hud (|1981li : these standard equations were implemented 
in the secular integrations to account for dissipative as 
well as non-dissipative tides. 

To illustrate the importance of non-secular effects for 
the dynamics of the systems under consideration, we 
show two examples in Figure [21 We consider a BH triple 
that is in the strongly non-secular regime defined by 
Equation m, as well as a stellar triple that is in the 
moderately non-secular regime defined by Equation (|9]). 
The evolution of the two systems was computed using 
both the secular equations of motion and ARCHAIN. 

The examples of Figure |2] demonstrate how the double- 
orbit-averaging procedure can lead to misleading results 
when applied to the dynamics of systems that evolve to 
attain very high orbital eccentricities. In the upper panel 
the three-body integration predicts that the inner binary 
will merge within a few LK oscillations whereas the secu¬ 
lar calculation predicts that the system does not merge in 
the considered time interval. More importantly, in the di¬ 
rect three-body integration the binary was found to have 
a substantial eccentricity of e, ~ 5 x 10“^ at the moment 
its peak gravitational wave frequency (Equation [13] be¬ 
low) reached 10 Hz. In the lower panel the stellar binary 
also attains higher eccentricities when integrated with 
ARCHAIN, which might induce a mass transfer event or 
reduce the circularization time due to tidal friction when 
compared to the secular integrations. 


3.2. Globular cluster models and initial conditions 


motion of arbitrarily tight binaries with arbitrarily large 
mass-ratio. The code also includes post-Newtonian (PN) 
non-dissipative IPN, 2PN and dissipative 2.5PN cor- 
rections to all pa i r-force s. We refer the reader to 
iMikkola fc Merriro (|2008ll for a more detailed descrip¬ 
tion of ARCHAIN. We also used the octupole level s ec- 
ular equations of motion given in iBlaes et al.l (120021) to 
evolve the systems, which allowed us to directly compare 
the predictions of the direct integrations to the orbit av¬ 
eraged results. 

In the case of triples containing an inner stellar bi¬ 
nary (stellar triples), to the PN terms we also added 


The initial conditions for the triple systems used in 
this analy sis were taken f rom t he Monte Carlo simu¬ 
lations of iMorscher et al.l (I20I5D . This study presents 
42 Monte Carlo dynamical simulations of massive star 
clusters containing large populations of stellar BHs. 
The Monte Carlo technique simulates two-body relax¬ 
ation in the Fokker-Planck approximation t hrough repre¬ 
sentative pair-wise sca ttering interactions ([Henonl [19711 : 
iChatteriee et al.ll2010D . 

Our imple mentation in the Cluster Monte Carlo 
(CMC) code (jPattabiraman et all I2013L and references 
therein) includes direct integration of strong three- and 
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four-body (binary-single and binary-binary) encounters, 
direct physical collisions, single and interacting binary 
stellar evolution, and tidal interactions with the Galaxy. 
Primordial triples are not included in these models. 
However, during strong binary-binary interactions, it 
is possible to for m stable hierarchical triple systems 
(|Rasio et al.lll99§l . Limitations in CMC currently re¬ 
quire that these triples be broken artificially at the end 
of the timestep. Nonetheless, whenever a stable triple is 
formed, its properties are logged, including the masses, 
radii, and star types for the stellar components, and 
the semi-major axes and eccentricities for the inner and 
outer orbits. Since we lack the information regarding 
the mutual orientation of the two orbits, we sample 
and cos(/) randomly from a uniform distribution with 
65° < I < 115°. The initial orbital phases in the three- 
body integrations were also randomly distributed. 

We extracted the properties for every stable triple sys¬ 
tem formed through dynamical interactions over 12 Gyr 
across all of the 42 models presented in Morscher et al. 
2015, and from five additional models provided by the 
same authors (Morscher, private communication). The 
m ain properties of the 4 7 cluster models are reported 
in iRodriguez et aH (j2015l l: these span a large range of 
masses (from 2 x 10® Mq to 1.6 x 10 ® Mq), of virial radii 
(from 0.5 pc to 4 pc), and metallicities (Z=0.0005, 0.0001 
and 0.005). We note that, because these triple systems 
are destroyed in the Monte Carlo simulation, it is possi¬ 
ble for the components of these triples to subsequently 
form new triple systems, when in reality they could in 
principle survive for a significant period of time. 

In the high-density stellar environment of GCs, triples 
may be perturbed through encounters with other pass¬ 
ing stars on timescales that can be shorter than the rel¬ 
evant LK timescale. Such encounters will alter the or¬ 
bital properties of the triple significantly, or even dis¬ 
rupt it. To account for this we set the final integration 
time in our simulations equal to the encounter timescale 
for collisions with o ther stars (|Binnev fc Tremainelll987l : 
llvanova et al.l[20?i^ . 


2 1 
Aftri + (AL) 


!8.5 X 10^^ yrPy 


1 -f 913 


zr'2 (JiQ 


( 12 ) 


where (Tiq is the central velocity dispersion of the cluster 
in units of 10 km s“^, Mtri = Mb + m 2 is the mass of the 
triple, (M) is the mass of an average star in the cluster, 
and rib is the local number density of stars in units of 
10® pc“®. The values of n, (M) and a in equation (fT^ 
were obtained directly from the Monte Carlo models and 
correspond to to their values in the cluster core at the 
moment the triple system was formed. We then integrate 
each triple until a time Tenc was reached, or until the 
inner binary had merged. 

In total, the Monte Carlo models contain 8864 triples. 
Of these 5238 are BH triples; 3626 have at least one non- 
BH component in the inner binary, and of these, in the 
inner binary, 2940 have two MS stars, 260 have a BH plus 
another MS star, 218 have a white dwarf plus another MS 
star, 4 have a neutron star (NS) plus a MS star, 163 have 
a MS star plus an evolved star (red giant or later type), 2 
have a BH plus an evolved star. The remaining systems 



Fig. 3. — Ratio of the triple survival time in the GC core to 
the LK timescale as a function of the ratio a 2 (1 — 62 )/^!. Tri angl e 
symbols represent systems that satisfy the condition Equation (I10[l . 
open circles are systems which satis fy the condition Equation 
but not the condition Equation Gnii and plus symbols are systems 
which satisfy neither of these two conditions. We evolved each of 
these systems ten times selecting randomly the initial inclination as 
detailed in the text. Blue colored points are systems in which in at 
least of the 10 realizations we found one merger within Tenc. Red 
points are systems in which for at least one of the ten realizations 
the inner binary eccentricity was larger than 0.9 at the moment its 
peak GW frequency became larger than lOHz. As predicted, most 
eccentric mergers occur in triples for which the orbit averaging over 
the inner binary breaks down. 

include 32 triples with an inner white dwarf binary, 5 
triples with an inner BH plus white dwarf binary, and 
2 triples with an inner BH plus NS binary. No triple 
with an inner NS binary is seen in any of these Monte 
Carlo models. Of the 5238 BH triples, 3852 satisfy the 
condition Tlk > Lenc, whereas 2519 of the 2940 stellar 
triples with an inner MS binary satisfied this condition. 


4. RESULTS: BLACK HOLE TRIPLES 

For each of the 3852 BH triples extracted from the 
Monte Carlo models we made ten realizations, each 
adopting a random orientation of the two orbits and ran¬ 
dom orbital phases. Of the total 38520 BH triples that 
we evolved with ARCHAIN, 2080 resulted in a merger 
of the two inner BHs. The total number of BH mergers 
obtained when using the orbit averaged equations of mo¬ 
tion was 1796 instead. We found that the merger time 
of such binaries after the triple was formed was typically 
very short and always shorter than < 10® yr. Thus, the 
BH mergers discussed in what follows will occur inside 
the GCs even when the recoil velocity at formation can 
be large enough to eject the triple from the cluster core. 
In this section we discuss the properties of the merg¬ 
ing compact binaries including their eccentricity, mass, 
merger time distribution, and event rates. 

Eccentric binaries emit a GW signal with a broad 
spectrum of frequencies, the peak gravitational wave 
frequency associated with the harmonic which leads to 
the maximal e mission of GW radiation can be estimated 
as (I WenI 1200311 


/gw = 


vCT/ (l + ei)ii9®4 


(13) 


r [ai(l-e?)]^'= 

We approximately follow the time evolution of the 
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Fig. 4.— Eccentricity distribution of merging BH binaries at the moment they first enter the lOHz (left panel) and 40Hz (right panel) 
frequency bands. While the octupole level secular equation of motions predict that only a few percent of systems will have finite eccentricity 
as they enter the aLIGO band, accurate body integrations show that 20% 10%) of BH mergers in GC will have an eccentricity 

larger than 0.1 at lOHz (40Hz) frequency. About 10% 5%) of all mergers will have extremely high eccentricities, i.e. 1 — e < 10“^^, at 

lOHz (40Hz) frequency. Note that the stippled regions are in front in both panels, which means that the lack of stippled regions at high 
eccentricities is because there are none, rather than because they are hidden behind the solid hystograms. 


GW frequency of the inner BH binaries through Equa¬ 
tion d. This allowed us to estimate the eccentricity 
of the BH binary when its GW frequency is > 10 Hz, 
i.e., when it wou ld be large enough t o be into the aLIGO 
frequency band (|Abadie et al.ll201^ . 

Figure [3] gives the ratio of the triple survival timescale, 
Tenc, to the LK timescale plotted against the value of 
02(1 — 62 ) 701 . From this plot we see that non-secular 
dynamical effects are expected to become important to 
the evolution of most BH triples in our models, possi¬ 
bly leading to the formation of eccentric GW sources. 
Also, most BH mergers are found to occur in moder¬ 
ately hierarchical triples with 02(1 — e 2 )/ai < 10. This 
is expected for at least two reasons: (1) the closer the 
outer body to the inner binary the less significant is the 
quenching of the LK cycl es due to relativis tic precession 
of the inner binary (e.g., iBlaes et aI1l2002D . so that typ¬ 
ically the smaller the ratio 02(1 — e 2 )/ai the larger the 
maximum eccentricity attained by the inner binary and, 
consequently, the shorter its merger time; (2) triples with 
small 02(1 — 62)701 ratio, have also larger Tenc 72 LK ra¬ 
tio, which naturally leads to a higher chance for a merger 
before the triple is disrupted by gravitational encounters 
with other stars. In Figure [3] we also identify those sys¬ 
tems for which at least one of the 10 random realiza¬ 
tions led to a merger (blue points) and those for which 
at least one realization produced a BH binary merger 
with an eccentricity larger than 0.9 at > 1(1 Hz fre¬ 
quency (red symbols). As predicted on the basis of our 
discussion in Section [2l most eccentric mergers occur at 
02(1 — 62 ) 7^1 ^ 10 , near the boundary for instability. 

Figure [1] shows the eccentricity distribution of all BH 
binary mergers in our simulations, when the associated 
/gw fi'i'st enters the 10 Hz and 40 Hz frequency bands. 
Figure [5] and [HI show the corresponding merger time and 
mass distribution of the merging binaries. The distri¬ 
butions showed in Figures IH [5] and [6] were obtained by 
weighing the number of mergers for each cluster model 
by the likelihood of that model to represent a typical GG 


in the Milky Way (MW). More in detail, the weights are 
obtained by creating a kernel density estimate (Wmw) 
of the MW GCs on the fundamental plane (which we 
take to be mass and ratio of the half to core radius), 
then estimating the weight for each model using the ker¬ 
nel density estimate at the position of that model on the 
fundamental plane. In this way, cluster models that are 
more likely to be drawn from the same distribution of 
MW GGs are more heavily weighted. The weight for a 
model of total mass Mqc , core radius Rc and half mass 
radius Rh after 12 Gyr of evolution is computed as: 


W{M, RjRh) 


Wmw(-^, Rc/Rh) 
WMode\s{M,Rc/Rh) 


(14) 


where we have divided by WModeis, the kernel density es¬ 
timate of the models themselves. This serves to normal¬ 
ize the distribution, so that regions of parameter space 
that are over -sampled by the models are given lower 
weights fsee iRodriguez et al.l 120151 for more details). 
The fraction of mergers with a given property (e.g., total 
mass, eccentricity) is then simply / = 
with Ni the number of mergers occurring in the Zth clus- 
t er model. _ 

iHuerta fc BrownI (|2013f) showed that for eccentricities 
less than 61 < 0.1 at « 10 Hz, circular templates will be 
effective at recovering the GW signal of eccentric sources. 
FigurelHshows that approximately 20% of all BH mergers 
in our three-body integrations had an eccentricity 6 > 0.1 
at > 10 Hz frequency. This percentage drops to ^ 10% 
at 40 Hz frequency. The difference with the results of 
the secular equations of motion is evident in these plots. 
The secular integrations clearly underestimate the num¬ 
ber of eccentric mergers producing just a few percent of 
inspirals with e > 0.1 at 10 Hz frequency. The direct 
three-body integrations also produce a significant pop¬ 
ulation (^ 15% of the total) of highly eccentric sources 
in the aLIGO frequency window, which are fully missed 
when evolving the triples with the secular equations of 
motion. These sources will start to inspiral due to GW 
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radiation energy loss within the aLIGO band with an 
extremely high eccentricity, 1 — ei ;< 10“^, and a charac¬ 
teristic GW frequency. In fact, the angular momentum 
of BH binaries that evolve into the strongly non secular 
regime can be arbitrarily small at the moment energy 
loss due to GW radiation starts to become dominant. 
For this reason the GW frequency at which the binary 
enters the aLIGO band, which depends on the binary 
angular momentum via Equation (1131) . could be signifi- 
cantly larger than 10 Hz. (see the example of Figure 3 in 
lAntonini et al.l[2014l l. We note that the detectability of 
such sources also depends on how much energy is radi¬ 
ated around 100 Hz where aLIGO is more sensitive. So 
it is important to see at what frequency these sources 
hrst enter the aLIGO sensitivity window. Based on our 
simulations we found that all eccentric BH binaries will 
enter the aLIGO band at ;< 300 Hz and 80% of them had 
/gw < 50 Hz when they f irst enter the aLIGO frequency 
band (see also Figure 7 in lAntonini et al.ll2014ll . We con¬ 
clude that most eccentric sources could be detected by 
aLIGO, provided that an efficient search strategy is used. 

Figure [5] displays the average merger time distribution 
of BH binaries induced by the LK mechanism in our mod¬ 
els. The figure shows a clear peak around IGyr of evolu¬ 
tion with a long tail extending to the present epoch. The 
declining merger rate is mostly due to the decreasing to¬ 
tal number of BHs in our models with time. As the rate 
depends on the total number of BHs in the cluster core 
it drops with time as BHs are continuously ejected from 
the cluster through strong dynamical interactions. 

Figure [S] displays the chirp mass, Mch = 
(mouii)^'^^ / (mo-l-as well as the total mass 
distribution of the merging BH binaries. These distribu¬ 
tions appear to be consistent with tha t of dynamically 
forme d binaries shown in Figure 2 of IRodriguez et all 
(|20I5H . Figure |6] shows a mass distribution peaked 
around Mch ~ 17 Mq, while the highest chirp mass of 
BH binaries formed by pure stellar evolution over all 
our models was ~ 13 Mq. This latter value for the 
chirp mass is the maximum for primordial binaries in 
our GG simulations that evolve d without sign i ficant 
dynamical encounters (see also iDowning et al.l 120101 : 
IRodriguez et al.l 1201511 . The fact that the typical 
chirp mass of the merging BH binaries in Figure [5] is 
substantially larger than r; 13 Mq is a result of the 
naturally stronger segregation of the heaviest BHs as 
well as the preferential ejection of the lightest BHs 
from the clusters during dynamical interactions. As the 
heavier BHs segregate more efficiently to the cluster 
center and the lighter BHs are ejected through binary 
single interactions, the former will tend to dominate the 
cluster core, where binary-binary interactions can lead 
to the formation of stable BH triples and to LK induced 
mergers. 

4.1. Black Hole Merger Detection Rate 

One of the advantages of using realistic cluster Monte 
Carlo models is that it allows us to make direct pre¬ 
dictions about the mass and merger time distribution 
of the coalescing BH binaries, which is fundamental to 
make reliable estimates for the event and detection rates. 
As we have complete information about the distribution 
of sources in time and chirp mass, we can compare our 



Merger Time (Gyr) 

Fig. 5.— Merger time distribution of BH binaries due to the 
LK mechanism in GCs. Continue lines give the merger time dis¬ 
tribution of all BH binary mergers, dashed line is for BH binaries 
that have e > 0.9 when /gw ^ 10 Hz. Insert panel gives the 
corresponding cumulative distributions. 



Chirp Mass, Total Mass (M^^^) 

Fig. 6. — Continue black lines show the chirp mass distribution of 
BH binary mergers due to the LK mechanism in CCs. Dot-dashed 
blue lines give the distribution of total mass for the merging bina¬ 
ries. Insert panel gives the corresponding cumulative distributions. 

ensemble of GC models to observations of MW and ex- 
tragalactic GCs and estimate the merger rate for a single 
aLIGO detector. 

To compute the rate of detectable sources per year 
from triple syst e ms, w e follow a similar procedure to 
IRodriguez et al.l (|2015ll . The rate is expressed as the 
following double integral over source chirp mass and red- 
shift: 

Rd = jj'Jl{Mc,z)fd{Mc,z)^^dMcdz. (15) 
where 

• TZ{Aic, z) is the rate of BH binary mergers with 
chirp mass M.c sX redshift z due to the LK mech¬ 
anism. 

• fd{-Mc, z) is the fraction of sources with chirp mass 
Ale at redshift z that are detectable by a single 
aLIGO detector. 

• dVrJd z is the co moving volume at a given red¬ 
shift (|Hogd [l9^ . We assume cosmological pa- 
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rameters of = 0.309, = 0.691, and 

h = 0.677, from the combined Planck results 
(IPlanck Collaboration et al.ll2015ll . 

• dts/dto = 1/(1 + z) is the time dilation between 
a clock measuring the merger rate at the source 
versus a clock on Earth. 

The rate is expressed as the product of the mean num¬ 
ber of triple sources per GC, the distribution of sources 
in Mc — z space, and the spatial density of GCs per Mpc^ 
in the local universe (0.77 Mpc“^, from iRodriguez et al.l 
(|2015f) . Supplemental Materials). Symbolically, this be¬ 
comes TZ{Mc,z) = (N) X P{Mc,z) X pGC- 

The components of 'lZ{Mc, z) are computed as follows. 
The mean number of sources per GG is found by per¬ 
forming a weighted linear regression between the num¬ 
ber of triple mergers a GC produces over its lifetime, 
and the final mass of the GC at 12 Gyr. The num¬ 
ber of triple mergers for a given cluster is obtained by 
multiplying the total number of mergers obtained from 
the three-body integrations by a factor cos(65°) x 0.1; 
where the first term takes into account the fact that 
we have only considered inclinations 65° < / < 115° 
and the second term the fact that for each triple we 
have performed ten integrations with random orbital ori¬ 
entations and phases. The weights are assigned using 
Wmw{M, Rc/Rh), the numerator from Equation (IT4ll . 
This linear regression is t hen multiplied by a universal 
GC luminosity function (jHarris et al.l 120141 1. and inte¬ 
grated over mass (assuming a mass-to-light ratio of 2) 
to yield a mean number of sources. This average is then 
multiplied by the distribution of sources, found by ran¬ 
domly drawing a number of mergers from each model 
according to the weights in Equation (|T4)l . then gener¬ 
ating a kernel density estimate of these sources in chirp 
mass/ redshift space. Note that unlike iRodriguez et al.l 
we do not separately consider the populations of 
low-metallicity and high-metallicity GCs. 

Einally, we model the efficiency of an aLIGO detector 
by calculating the fraction of sources at chirp mass A4c 
and redshift z that aLIGO could detect, marginalized 
over all possible sky locations and bi nary orientations. 
We use IMRPhenomC waveforms (iSantamarfa et aO 
[Mnli and the projected zero-detuning, high-power 
aLIGO noise curve 0. Note that this template fam¬ 
ily does not consider the eccentricity of sources, which 
could dramatically alt er the gravitational waveforms 
(|Hnerta fc BrownI 120131 1. Such simplification does not 
significantly alter the total detection rate estimates pre¬ 
sented here; in fact, our simulations show that at 10 Hz 
frequency ^ 80% of coalescing BH binaries have an ec¬ 
centricity ei ^ 0.1, which is small enough that circular 
templates will be efficient at recovering their GW si gnal 
(jBrown fc Zimmermanl[2?)l(Tl : iHuerta fc Brownll2013[l . 

As P{A4c, z) is dependent on the random draw of 
mergers, we report the mean of Rd over 100 separate 
draws. We also provide optimistic and pessimistic rate 
estimates to approximate the uncertainties associated 
with our assumptions. The optimistic rates are com¬ 
puted by assuming the +la value of the 100 draws of 
Rd, the +la uncertainty for the linear regression between 

^ LIGO Document T0900288-v3 


GC mass and number of mergers, a highly optimistic up¬ 
per limit on the GC mass function (2 x 10^Mq), and the 
upper -limit on pcc of 2.3 Mpc~^ from IRodriguez et al.l 
(1201511 . Conversely, the pessimistic rate estimate as¬ 
sumes the corresponding — la uncertainties, a pessimistic 
upper limit on the GC mass function (4 x 10^Mq), and 
the lower-limit of pcc of 0.32 Mpc“^. 

Eor all LK driven events, we find that a single aLIGO 
detector could detect ~ 0.5 events per year (pessimisti¬ 
cally ^ 0.1 yr“^, optimistically ^ 2 yr“^). 

We performed the same analysis described above for 
all eccentric sources produced in our three-body simu¬ 
lations, i.e., we computed the detection rate solely for 
eccentric systems assuming perfect templates for detect¬ 
ing eccentric binaries. Thus, the results of this computa¬ 
tion do not represent the number of sources that aLIGO 
could detect with current searches. Instead, they suggest 
an estimate of the number of sources that aLIGO could 
detect if optimal matched-filtering searches for eccentric 
sources were to be used. For events with eccentricities 
greater than 0.1, this rate is ~ 0.2 yr“^ (pessimistically 
^ 0.02 yr“^, optimistically ^ 0.5 yr“^), while for merg¬ 
ers with eccentricities greater than 0.9, the rate becomes 
~ 0.2 yr“^ (pessimistically ~ 0.01 yr“^, optimistically 
~ 0.4 yr“^). We note also that that the correspond¬ 
ing waveforms of most eccentric GW sources discussed 
here are characterized by a sequence of bursts of energy 
emitted near periapsis rather than a continuous wave¬ 
form. Because of this, matched hltering techniques could 
be unpractical and a power stacking algorithm might be 
employed instead (jEast et al.l 1201311 . The disadvantage 
is that this will not give an optimal signal-to-noise ra¬ 
tio which might somewhat reduce the detection rates for 
eccentric inspirals we find. 

Table 1 summarizes the main results of our rate com¬ 
putation giving the predicted aLIGO detection rate of 
binary BH mergers from GCs. This table also gives the 
corresponding compact binary coalescence rates per GC 
per Myr. Table 1 gives the rate of LK induced BH merg¬ 
ers (hierarchical triples), and the rate of mergers due to 
the LK mechanism that have an eccentricity larger than 
0.1 at 10 Hz frequency (eccentric inspirals). From our 
Monte Carlo models we also computed the merger rate of 
hard escapers that will merge in the galactic field within 
one Hubble time (ejected binaries) and the rate of merg¬ 
ers occurring inside the clusters that are not due to the 
LK mechanism (in-cluster mergers). These latter include 
mergers due to binary-single dynamical interactions, as 
well as those mergers due to binary stellar evolut ion and 
single-single captures (see IRodriguez et al.l[2015l for de¬ 
tails). 

In agreement with previous calcu lations (e.g., 
lO’Learv et al .11200^ iDowning et aLll^OlOH we find that 
the total event rate is largely dominated by outside 
cluster mergers. We stress again that the BH mergers 
due to the LK mechanism will occur inside the GC. 
Thus, by comparing the rate estimates for triples to the 
rate of in-cluster mergers due to other processes in Table 
1 we can conclude that the detection rate of BH mergers 
occurring inside the cluster is dominated by the rate 
of LK induced mergers. This is somewhat surprising 
given that in-cluster mergers that do not occur in triples 
account for approximately 10% of the total binary BH 
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TABLE 1 

aLIGO detection and merger rates of coalescing BH binaries from GCs. (In-gluster mergers are defined here as all events occurring inside 
THE GG Monte Carlo models that are not due to the LK mechanism. As described in the text, the detection rate for eccentric sources 
( ei[> 10 Hz] > 0.1) REPRESENTS THE NUMBER OF INSPIRALS THAT COULD BE DETECTED WITH DEDICATED SEARCH STRATEGIES.) 


Source 

low 

Detection rates (yr 
realistic 

high 

low 

Merger rates (GC ^ Myr 
realistic 

high 

Hierarchical triples 

0.11 

0.47 

1.93 

1.9 X lO-'^ 

2.8 X lO-'^ 

3.3 X lO-'^ 

Eccentric inspirals 

0.021 

0.20 

0.48 

3.4 X 10“® 

5.7 X 10“® 

7.3 X 10"® 

In-cluster mergers 

0.011 

0.062 

0.60 

1.9 X lO-'^ 

2.8 X 10“® 

3.4 X 10“® 

Ejected binaries 

10 

30 

100 

0.015 

0.023 

0.027 



1 10 100 1000 

02(1-62)70^ 

Fig. 7. — Similar to Figure|3]but for stellar triples. Blue symbols 
are systems In which one component of the inner binary crossed 
the Roche limit during the evolution. 

mergers over 12Gyr. The low aLIGO event rate we 
find is a consequence of the fact that the majority of 
these mergers occur early in the cluster’s lifetime, and 
are relatively low-mass. Since aLIGO is less sensitive 
to low-mass, high-redshift sources, the rate of detected 
sources drops. We conclude that mergers induced by 
LK oscillations in hierarchical triples are likely to be 
the main contributor to in-cluster gravitational wave 
sources for aLIGO detectors, and represent ^ 1% of the 
anticipated total detection rate of BH binary mergers 
assembled in GCs. 

5. RESULTS: STELLAR TRIPLES 

In total, our Monte Carlo models produced 2940 triples 
in which both components of the inner binary were MS 
stars (i.e., stellar triples). Following the same proce¬ 
dure adopted above, we made ten realizations for each 
triple, taking a random orientation between the two or¬ 
bits and random orbital phases. We then integrated the 
triples forward in time using ARCHAIN if the condition 
([HI) (with Udiss = 2[i?o + ^ 1 ]) was satisfied; the rest of 
the initial conditions were integrated using the octupole 
level secular equations of mot ion based on the dou¬ 
ble average Hamiltonian fe.g.. Il^oz fc: Fabrvckvl [2014 
iMichaelv fc Pere't^l20I4l) . Each integration was termi¬ 
nated either at Tenc, or when one of the binary compo¬ 
nents had passed within its Roche limit. 

In many of our simulations the inner binary reaches 
very high eccentricities, implying a high probability that 
the stars will cross each other’s Roche limit and transfer 
mass. To account for the possibility of crossing the Roche 


limit, we check whether the condition dijLn^ij < Ri (or 
ai(I — < Ri in orbit average integrations) was 

met, with dij the instantaneou s distance betwe en the two 
inner binary components and (|EggletonllI983ll : 


0.6 -b In ^I -I- 

( 16 ) 

As often done in the literature, we stop the integration 
and assumed that the binary had merged if one of the 
two b inary components had p assed within its Roche limit 
le.g.. iNaoz fc Fabrvckvl [20I4I1 . We found that in 10327 
out of the overall 29400 triples, one of the inner binary 
components crossed its Roche limit. After accounting for 
the fact that we have explored a limited range of orbital 
inclinations, we find an overall probability of ~ 10% that 
a stellar triple assembled dynamically in a GC will lead 
to a mass transfer event and possibly to a stellar merger. 

Figure [7] gives the ratio of the triple survival timescale, 
to the LK timescale plotted as a function of 02(1 — 
e 2 )/ai. Mass-transfer systems are indicated by blue col¬ 
ored symbols. This plot shows that the survival time- 
scale for these systems can be many orders of magnitude 
larger than the corresponding LK time-scale. The ratio 
Tenc/^LK for these systems is also typically larger than 
that for BH triples (compare to Figure [3|), reflecting the 
different tertiary period distribution of the two popula¬ 
tions shown in Figure[T]and the larger binary mass in the 
case of BH components. In Figure[7]we identify those sys¬ 
tems for which the orbit average approximation breaks 
down. In agreement with our previous discussion, most 
of these systems correspond to compact configurations 
in which the tertiary comes closer to the inner binary 
than a separation 02(1 — e 2 )/ai < 100. We also identify 
those systems in which in at least one of the ten random 
realizations one of the inner binary components crossed 
its Roche limit. Not surprisingly, close outer orbits are 
more likely to result in a mass transfer event between the 
inner binary components. 

In Figure [8] we consider the distribution of orbital pe¬ 
riods. Figures 0 and Figure [10] display respectively the 
distribution of merger times and mass of the merger rem¬ 
nant for systems that crossed their Roche limit. These 
distributions were obtained using the weighing technique 
described in Sectional i.e., weighing the number of sys¬ 
tems for each cluster model by the likelihood of that 
model to represent a MW-like GG. 

In Figure |5| we show the initial and final orbital period 
distributions of the systems that did not cross the Roche 
limit, as well as the initial period distributions of the sys¬ 
tems that crossed the Roche limit during the evolution. 
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Fig. 8. — Initial (black lines) and final (dot-dashed red lines) period distributions for systems in which the inner binary components 
did not cross the Roche limit during the evolution. Dashed blue lines show the initial period distribution of systems that crossed their 
Roche limit during the evolution. The distributions are normalized to the total number of systems that were evolved, such that the sum 
of the integrals of the distribution of not merging systems plus the distribution of merging systems is unity. As explained in the text these 
distributions should be interpreted as typical distributions for MW-like GCs. 


As shown in this figure the peak of the period distribu¬ 
tion of the surviving binaries shifts from ~ 3.2 days to 
^ 1 day. This is a consequence of energy loss due to tidal 
dissipation that the binary experiences during the high 
eccentricity phases of LK cycles. If during the evolution 
the timescale associated with precession due to tidal and 
PN terms becomes shorter than the LK timescale, the ec¬ 
centricity oscillations due to the interaction with the ter¬ 
tiary are suppressed. In this case, after many LK cycles 
tidal dissipation can shrink the binary orbital semi-major 
axis by a large factor, causing the binary to effectively 
decouple from the tertiary. 

Figure m shows that although merging systems are pro¬ 
duced at any value of Pi and P 2 , they are more likely to 
originate from systems with the tightest outer orbits and 
widest inner binaries. In fact, about half of the systems 
with P 2 < 100 days or with Pi > 10 days led to stel¬ 
lar mergers. Generally in order for the inner binary to 
attain a quasi-radial orbit such that its components can 
merge, its eccentricity has to be excited on a timescale 
much shorter than the typical extra precession timescale 
(such as tides and relativistic precession). It can be easily 
shown that the effectiveness of such extra so urces of pre¬ 
cessio n decreases with the ratio 02/01 (e.g., IBlaes et all 
I 2 OO 2 D , which is consistent with the fact that mergers are 
more likely to occur for wide inner binaries and close 
outer orbits. 

The distribution of merger times for stellar binaries is 
shown in Figure |9l Given that the merger time through 
LK cycles is generally very short, the number of merging 
systems at a given time is also a good measure of the 
number of hierarchical triples that are formed dynami¬ 
cally in the cluster core at that time interval. Figure IHl 
shows that ~ 40% of all stellar mergers occur in the first 
Gyr of cluster evolution. This is expected given that it 
is only during these early times that the stars dominate 
the cluster core where stellar binary-binary encounters 
can take place at a high rate and lead to the assembly 
of stellar triples. After ^ I Gyr of evolution, the BHs 
have already decayed to the center through dynamical 
friction becoming the dominant population in the clus¬ 


ter core. As the BHs become dominant, interactions in¬ 
volving stellar binaries become less frequent, naturally 
reducing the formation rate of triples with stellar com¬ 
ponents and favoring the formation of BH triple systems 
instead. 

Finally, Figure ITUl displays the total mass distribution 
of all binaries that underwent Roche-lobe overflow dur¬ 
ing their evolution. This plot shows that the merging 
systems have masses up to ~ 8 Mq with a peak at about 
I Mq. The most massive mergers, including a few at 
> 4 Mq , were found to occur only very early in the first 
few Gyrs. As illustrated by the dashed-blue lines in Fig¬ 
ure [ini no merger with mass larger than ^ 2 Mq was 
found during the last 4 Gyr of cluster evolution. These 
late mergers have a typical mass of 0.5 — I Mq, as ex¬ 
pected being about twice the mass of an average MS star 
of an old single stellar population. 

5.1. Formation of blue stragglers 

BSSs are main sequence stars that appear to be hot¬ 
ter and more luminous than the turn-off point for their 
parent cluster population. As such they appear to be 
younger than the rest of the stellar population (e.g., 

I'LSellT^ . 

BSSs are created when a MS star is replenished with 
new supply of Hydrogen in its core. Hence, they can re¬ 
main in the MS longer than expected of an undisturbed 
star. This can happen through mass transfer in a binary 
syst em or via m erger/collision with a t least one MS star 
seelPeretsll20I5l for a recent review). iPerets fc Fabryck"^ 
200911 suggested that triples might be a natural progen¬ 
itor for BSSs, as the high eccentricities attained dur¬ 
ing LK cycles can lead to tidal interaction, mass loss 
or even a merger of the inner binary components. Al¬ 
though dynamically formed triples are expected to form 
quite efficiently in GCs, their long term dynamical evo¬ 
lution in GCs has not been yet investigated. Hence, until 
now their role in BSS formation remains largely uncon¬ 
strained. In what follows we use our Monte Carlo models 
coupled with the results of the three-body integrations 
to determine the relative contribution of triples to BSSs 
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in GCs. 

Besides BSSs possibly produced by the LK mechanism, 
we also computed the number of BSSs due to binary- 
mediated interactions, i.e., mass transfer and physical 
collisions. We computed the number of BSSs at the ob¬ 
servation time of ~ 12 Gyr in each model in the follow¬ 
ing way. We first estimate the MS turn-off mass (ttito) 
given the metallicity and age of the model cluster in ques¬ 
tion. From the snapshot containing all stars we then 
extract those that are still on their MS and has mass 
> 1.05 X ttito- In theory, any MS star with mass above 
m-TO are BSSs, however, observationally, it is hard to 
distinguish BSSs from the MS if their masses are not 
high enough compared to mxo ■ This simple mass-based 
prescription results in very good agreement between ob¬ 
servationally extracted BSS numbers, their stellar prop¬ 
erties, and known correlations betwee n BSS and observ¬ 
able properties of the cluster (e.g.. [Leigh et al.l 120111 : 
iChatteriee et al.ir2013l : iSills et al.ll2013t for more details'). 

Figure [TT] shows the relative importance of BSS pro¬ 
duction via traditional binary-mediated channels and 
triples. Since the triples were not followed in the Monte 
Carlo models, BSSs created via LK cycles leading to mass 
transfer or merger between the inner binary can only be 
calculated for each triple produced in the models. For 
these cases we cannot directly consider the subsequent 
stellar evolution of the BSSs. For example, whether fol¬ 
lowed by a merger the BSS created via a triple-mediated 
channel later evolves off of the MS before the observa¬ 
tion time of ~ 12 Gyr, cannot be directly determined 
in this setup. Indeed, in detailed models it had been 
found that BSSs can live in the MS for several billion 
years after formation. It is also expected that the me¬ 
dian age since formation for the observed BSSs can vary 
depending on the domina nt formation channel f or the 
BSSs in a real cluster fe.g.. IChatteriee et al.ll2013l) . The 
real age since formation for BSSs, even in these detailed 
models might not be accurate, since the residual reju¬ 
venated lifetime for a real BSS depends intricately on 
the details of the degree of H-mixing in the core which 
in turn depends on the details of the interaction that 
produced it (e.g., ISills et al.ll2001l:lLombardi et al. 1120021: 

iChen &: Han ll2009HChatteriee et al.ll2013[l . Instead, we 
report two quantities for the triple-mediated BSSs- 1) 
The total number of mergers or Roche-filling architec¬ 
tures created from dynamically formed triples during the 
full evolution of the model clusters, and 2) the same, 
but for binaries with a total mass > 1.05 x mxo and 
created only within the last 4 Gyr, an estimated upper 
limit of today’s observ ed BSSs in typical Galactic GCs 
(IChatteriee et al.]l20I3ll . We hnd that the relative im¬ 
portance of triple-mediated channels is low compared to 
binary-mediated channels for BSS production. The im¬ 
portance of triple-mediated channels is < 10% for the 
cluster models considered in this work. 

There is a lot of interest in BSSs found in binaries. 
Often these systems are thought to be produced via 
stable mass transfer in a binary where the donor stays 
bound to the accretor. On the other hand, these may 
also be produced via triple-mediated channels where 
the inner binary merges a nd the outer binary is hard 
and remains unbroken (e.g . dFabrvckv fc Tremainell2007l: 
Perets &: Fabrvek^ 120091 : iPeretsI 1201^ iGosnell et al l 
20I5|1 . We find that the contribution from triple- 



Merger Time (Gyr) 

Fig. 9. — Merger time distribution of stellar binaries produced 
by the LK mechanism in a MW-like GC. Insert panel shows the 
corresponding cumulative distribution. The distribution has been 
normalized to the total number of merging systems. 



Total Mass 

Fig. 10. — Distribution of the total mass of binaries that crossed 
the Roche limit during the evolution. Assuming that these bina¬ 
ries will end up merging, such distributions can be interpreted as 
those corresponding to stellar merger products induced by the LK 
mechanism in a MW-like GC. The dashed-blue lines show the mass 
distribution of systems that merged in the last 4Gyr of the cluster 
evolution. These distributions have been normalized to the total 
number of merging systems. 

mediated channels for binary BSSs is also small being 
< 10% of the number of mergers due to binary-mediated 
processes. In addition, we find that the binary BSSs com¬ 
ing via binary-mediated channels are not simply coming 
from mass transfer in a binary. Collisionally produced 
BSSs can frequently acquire a binary companion via 
binary-single scattering encounters. Since BSSs are more 
massive compared to typical cluster stars at late ages, ini¬ 
tially single BSSs, created earlier via a merger or physical 
collision can frequently acquire a binary companion via 
exchange with another normal star in a binary through 
binary-single interactions. This is consistent with earlier 
work that investigated the detailed history and branch¬ 
ing ratios of various binary-mediated fo rmation channels 
for B SSs produced in similar models (|Ghatteriee et al.l 

6 . DISGUSSION 

6 .1. Implications for GW detectors and Comparison to 
other sources of eccentric GW inspirals 
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Fig. 11.— Contribution of triples to the population of BSSs (left panels) and binary BSSs (right panels) in GCs. In the upper left panels, 
filled and open circles represent respectively the total number of BSS candidates and stellar mergers formed in triples. The filled and 
open triangles show respectively the total number of BSSs and stellar mergers in our cluster models that were not produced in hierarchical 
triples. In the right panels filled circles represent BSSs from hard triples with MS star tertiary. These latter are the systems that will 
survive as binaries up to the present epoch and could be observed as binary BSSs. In the lower panels the filled circles give the ratio of the 
number of (binary) BSSs formed in triples (Ntr) to the number of BSSs due to other channels (Noe)- The open circles give the number of 
stellar mergers from triples devided by the total number of stellar mergers from binary mediated channels. These quantities are all plotted 
against the central density of the model at 12 Gyr. Note that a few clusters are in a state of core-collapse at this time, being characterized 
by very high central densities (pc > lO®M 0 pc“®). 


The first GW signal from a compact object binary 
is likely to be detected in the coming yea rs by Ad¬ 
vanced ground-based laser interferometers (|Aasi et al.l 
120151: lAcernese et al.l 1201 51) . Most GW searches adopt 
matched filtering techniques in which the detector sig¬ 
nal is cross-correlated with banks of theoretical gravita¬ 
tional waveforms to enable detection of the weak signals. 
Mainly due to computational limitation, such banks of 
theoretical waveforms are exclusively composed of cir¬ 
cular binaries. However, there are many reasons to 
extend the searches to also include eccentric inspirals. 
First, there have been recent developments of dedicated 
techniques that wi ll potentially en a ble the detection of 
eccentric sourc es (iTai et al.l 12014 iHuerta et al.l 12014 
iGoughlin et al.l[2015HAkcav et al.ll2015D . Second. eccen¬ 
tric inspirals produce a GW signal which is distinct from 
that of circularly inspiraling binaries, and provides more 
insights on th e strong-held dynam ics (jLoutrel et al.ll20l4 
iLe Tiedl2015t l^cav et all 1201511 . Finally, as discussed 
below, there are several dynamical mechanisms that can 
lead to the formation of eccentric GW sources. The de¬ 
tection of such eccentric sources could provide with im¬ 
portant information about properties (e.g., triple/binary 
fraction, central densities) of compact object populations 
in GGs and galactic nuclei that will be otherwise inac¬ 
cessible. 

Several mechanisms have been proposed for the for¬ 
mation of eccentrically inspiraling compact binaries for 
high frequency gravitational wave detectors. These in¬ 
clude: (i) BH-BH scattering in steep density cusps 
around massive black holes (MBHs: lO’Learv et al.l[200^ 
iKocsis &: Levinl[2012D . (ii) bi nary-single coinpact o bject 
encounters in star clusters (|Samsing et al.l 1^014) and 


(iii) LK resonance in hierarchical triple systems form- 
i ng in the dense stellar environment of galactic nuclei 
(lAntonini fc PeretsI l2012t) or stellar clusters (e.g.. iWei] 
120031 this paDerb In this section we discuss the rate for 
eccentric compact binary coalescences which could be de¬ 
tectable by aLIGO. With the caution that rate estimates 
remain subject to significant uncertainties, we find that 
BH binary mergers from hierarchical triples in GCs are 
likely to dominate the production of eccentric inspirals 
that could be detectable by Advanced ground-based laser 
interferometers. 

Single-single captures. In galactic nuclei with 
MBHs, nuclear relaxation times can be often less than 
a Hubble time, and can result in the formation of steep 
density cusps of stellar-mass BHs as a consequence of 
mass segregation against th e lower mass stars (e.g., 
iHonman fc Alexandeill200^ . lO’Learv et al.l (j2009f ) pro¬ 
posed that in such dense population environments BH bi¬ 
naries can efficiently form out of gravitational wave emis¬ 
sion during BH-BH encounters. After forming, these BH 
binaries rapidly insp iral and merge within a few hours. 
lO’Learv et al.l (I2009D found that ^ 90 % of these co¬ 
alescing binaries will have an eccentricity larger than 
0.9 inside the aLIGO frequency band. The predicted 
rate for eccentric c apture of IOM 01 BH-BH encounters is 
~ 0.01vr~^G pc~^ (lTsang|l2013D. 

Table 1 of lO’Learv et al.l '( 20091) presents the merger 
rate per MW-like galaxy and for different cusp models. 
The plausible pessimistic, likely, and optimistic rates can 
be taken from models A/33, E-2, and F-1 of their Ta¬ 
ble 1. The most optimistic scenario, given by model F-1, 
produces a merger rate of 1.5 x 10“^ Myr“^. The re¬ 
alistic scenario, identified here with model E-2, predicts 
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1.3 X 10“^ Myr“^. The pessimistic scenario, given by 
model A/33, produces 2 x 10“^ Myr“^. We can compare 
these merger rates with those corresponding to eccentric 
inspirals from GCs. After multiplying the merger rates in 
our Table 1 by a factor 200 (ro ughly the number of GCs 
in the MW; Harris et al.ll201 4fl. we see that the realistic 
merger rate for 10 ’Leary et al.l (|2009ll scenario appears to 
be about one order of magnitude smaller than the real- 
is tic merger ra t e of e ccentric inspirals from GCs. 

iHong fc Led (j2015l l carried out A^-body simulations of 
star clusters around MBHs to study the formation of 
BH binaries in nuclear clusters and to put constraints 
on th e aLlG O detection rates for the 10 ’Leary et al.l 
(|2009ll model. iHong fc Led (I2015D derived an aLlGO ex¬ 
pected detection rate of 0.02 — 14 yr“^ depending on 
the maximum horizon distance assumed and the mass 
rat io of MBH t o the s urrounding cluster. As also note 
by IHong fc Leel l)2015fl , these rates are likely to be an 
overestimate of the real event rate from BH-BH cap¬ 
tures. They were in fact obtained by assuming that the 
mass fraction of BHs is 44% of the total cluster mass, 
which is only reasonable within about one-tenth of the 
MB H influence radius of a mass -segregated nuclear clus¬ 
ter (|HoDman fc Alexanded[200^ . Given that most of the 
captures were found to occur near the half-mass radius 
of the nuclear cluster, where mass segregation will not 
significantly a lter the initial n umber fraction of different 
mass species (IAntoninill2014f) . a more reasonable choice 
would be to simply adopt a BH mass fraction similar to 
that expected on the basis of the initial mass function of 
the underlining population. For example, in the Galactic 
center there is evidence for a standard Kroupa-like initial 
mass function (jLdckmann et al.ll20ir)fl which will result in 
1% o nly of the total mass in BHs ( Hopman fc Alexander! 
1200611 . Using this standard choice for the initial mass 
function and adoptin g the same quadrati c scaling on the 
BH mass fraction of IHong fc Led (j2015[ l will cause the 


estimated detection rates to drop by approximately 3 or¬ 
ders of magnitude. 

Binary-single encounters. The in-cluster detection 
rate reported in Table 1 includes BH mergers due to 
single-binary interactions. Assuming that 1% of these 
merge rs will have a finite e ccentricity within the aLIGO 
band (jSamsing et al.l[20T^ results in a negligible contri¬ 
bution of binary-single encounters to the the population 
of eccentric inspirals when compared to the triple chan- 
neh _ 

iSamsing et al.l (|2014ll studied binary-single stellar scat¬ 
tering occurring in dense star clusters as a source of ec¬ 
centric NS-NS inspirals for aLIGO. During the chaotic 
gravitational interaction between the three bodies, a pair 
of compact objects can be driven to very high eccentric¬ 
ities such that they can i nspiral through gravi tational 
wave radiation and merge. ISamsing et al.l (l2014f) argued 
that 1% of dynamically assembled NS-NS merging bi¬ 
naries will have a substantial eccentricity at > 10 Hz 
frequency. ISamsing et al.l (j2014ll also give a simple order 
of magnitude estimate of the merger rate for dynami¬ 
cal NS-NS star eccentric inspirals assembled in GCs of 
0.7yr-iGpc-3. 

To compare to the event rate we found for eccen- 
tric LK induced m ergers in GCs we recompute the 
ISamsing et al.l (|2014l) rates by rescaling with the different 


GC nu mber density used in Section H?T] ISamsing et al.l 
(|2014f) calculate the rate of eccentric binaries in the 
aLIGO band adopting poc = 10 GC Mpc“^. Taking 
Pgc = 0.77 GC Mpc“^ instead, assuming that aLIGO 
will see NS- NS mergers up to a sky-averaged distance 
of 400 Mpc (|Abadie et al.l]2H^ . and optimistically as¬ 
suming that mergers are distributed uniformly over the 
lifetime of the globulars (as opposed to happening early 
on in the lifetime of the cluster), this translates into a 
possible aLIGO detection rate for NS-NS eccentric inspi¬ 
rals of ~ 0.02yr“^. 

We note that iCoughlin et al.l (|2015f) give in their Ta- 
ble H the potential detection rates for the model of 
ISamsing et al.1 (j2014l l. These rates appear to be larger 
by one order of magnitude than our rate estimate 
ab ove. We believe how ever that the rates computed 
bv ICoughlin et ahl (|2015f) were mistakenly enha nced due 
to a misinterp r etatio n of ISamsing et al.l (|2014fl results. 
ISamsing et al.1 (j2014ll found that about “1% of dynam¬ 
ically assembled non-eccen tric binaries” w i ll hav e a fi¬ 
nite eccentricity at 10 Hz. ICoughlin et al.1 (j2015tl com¬ 
puted the rates in their Table H as the 1 % of the total 
compa ct merger rate given in Table V of lAbadie et al.1 
(|2010H which corresponds, however, to predictions for 
field mergers. In order to co rrect for this we have to 
multiply the rates in Table H of iGoughlin et al.l (j2015f l by 
the expected fraction of all compact object mergers that 
are dynamically assembled in star cluste rs. Taking this 
fraction to be 10% of all NS -NS mergers (|Grindlav et al.1 
120061 : ISamsing et al.1 [20Tll . and assuming that 1% of 
these are eccentric inspirals results in 4 x 10“^yr~^, 
0.04yr“^, and 0.4yr“^ for the low, realistic, and high 
detection rates for binary NSs. These estimates appear 
to be in good agreement with the rate of ^ 0.02yr“^ we 
previously computed. 

Although very simplihed, the simple calculations pre¬ 
sented above show that, even when assuming a uniform 
rate of mergers, the rate of eccentric NS binaries in the 
aLIGO frequency band is likely 1-2 orders of magnitude 
smaller than the realistic estimates for eccentric BH bi¬ 
naries forming through the LK mechanism in GGs. 

Massive black hole me diated compact-object 
mergers in galactic nuclei. lAntonini fc Pereti (|2012f ) 
studied the evolution of compact object binaries orbiting 
a MBH. Near the galactic center, where orbits are nearly 
Keplerian, binaries form a triple system with the MBH 
in which this latter represents the outer perturber of the 
triple. The MBH can subsequently drive the inner binary 
orbit into high eccentricities through LK resonance, at 
which point the compact objec t binary can efficiently co - 
alesce through GW emission. lAntonini fc Pereti (|2012ll 
found that ~ 10% of such coalescing binaries will evolve 
through the dynamical non-secular evolution similar to 
that described in Section O leading to finite eccentricity 
s ources in the aLIGO sensit ivity window. 

lAntonini fc Pereti (|2012ll estimated a BH-BH ec¬ 
centric inspiral rate for a MW-like galaxy between 
10“® Myr“^ and 10“^ Myr“^. These rates are some¬ 
what comparable to the ra tes for BH-BH merge rs due 
to GW captures derived bv lO’Learv et all (j2009fl : how¬ 
ever, they are still smaller than the rates for eccentric 
inspirals assembled in GCs for a MW-like galaxy given 
in Table 1, and, even under the most optimistic assump- 
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tions, they are still only 3% of the realisti c and 0.01% of 
the high/optimistic rate in lAbadie et al.l (|2010r) . 

Given the re alistic BH-BH detec tion rate with aLIGO 
of ~ 20yr“^ (jAbadie et alJ l20i(ill . the predicted over¬ 
all rate of eccentric compact binary coalescence in the 
aLIGO band induced by M BH mediated LK cyc les in 
galactic nuclei is ~ 0.02yr“^ (|Coughlin et al.l[2015ll . We 
conclude that, the overall contribution of eccentric GW 
sources from galactic nuclei is likely negligible compared 
to the rate of eccentric mergers arising from hierarchical 
triples in GGs. 

6.2. Implications for blue straggler formation 

The BSSs observed abundantly in all clusters have long 
been identified as a way of identifying important stages 
in the host cluster’s history. However, interpretation of 
observed trends between the BSS number and cluster 
properties is dependent on the relative importance of the 
various BSS formation channels in these clusters. Espe¬ 
cially, there has been a long history of carefully consid¬ 
ering correlations between the number of BSSs produced 
in a cluster and several dynamically important observed 
cluster properties, including the total mass, total number 
of binaries, and th e central collision/interaction rate, and 
what they inean (iPiotto et al.l [200^ iLeigh et al.l 120071 : 

iKnigge et al.ll2009HChatteriee et al.ll2013|l . Although, it 
has been pointed out that formation of BSSs from triple- 
mediated iritera£Bonsmn^^e important, at least in open 
clusters (|Leigh fc Silld I201lll . this has never been self- 
consiste ntly investigj ated in realistic and evolving cluster 
models (|Peretsll2015ll . As part of this study we have con¬ 
sidered evolution of triples created dynamically in the 
cluster models and analyzed them as potential BSS pro¬ 
genitors. Upon comparison of the estimated numbers of 
the BSSs resulting from triple-mediated formation chan¬ 
nels to those formed via binary-mediated channels (mass 
transfer and physical collisions due to binary-mediated 
interactions), we find that, at least for GG-like densities, 
contribution from triple-mediated channels is quite low. 

For the triple-mediated BSSs, we did not consider 
whether upon formation the rejuvenated MS star will re¬ 
main in the region expected of BSSs in a color-magnitude 
diagram. The rejuvenated lifetime of BSSs is model de¬ 
pendent and hard to estimate. For the binary-mediated 
BSSs we assume full-mixing of Hydrogen. In case of 
triple-mediated BSSs, we simply compare the total num¬ 
ber of BSSs created throughout the lifetime of the GC 
and that formed within the last 4 Gyr and with mass 
> 1.05 X mxo- Thus, the actual contribution from triple- 
mediated interaction may be even lower if many of these 
BSSs evolve off of their rejuvenated MS before the ob¬ 
servation time of ^ 12 Gyr. 

We note that above we have not considered BSS for¬ 
mation through binary evolution in short period binaries 
affected by the LK mechanism coupled with tidal friction. 
The shrinkage of the inner binary semi-major axis as well 
as the high eccentricities induced via the LK mechanism, 
combined with an expansion in the size of each inner bi¬ 
nary star, could lead to coalescing/strongly-interacting 
post-MS binaries. In order to address the importance of 
this alternative channel we determined whether the stars 
would undergo mass transfer by requiring their periapsis 
distance to be smaller than the Roche limit computed 
through Equation (fT6l) and setting the stellar radius of 


the most massive component e qual to that evalua ted at 
the tip of the red giant branch ijHurlev et al.ll2000l) . The 
periapsis distance was evaluated using the orbital param¬ 
eters at the end of the three-body integrations. We only 
considered stars with mass >0.7 Mq so that their MS 
life time was shorter than the host GC age. Also, we 
did not consider systems that will have undergone mass 
transfer on the red giant branch even without the help of 
LK cycles and tidal dissipation. We found that only 694 
of the total 19073 binaries that did not merge during the 
MS would merge on the red giant branch through a com¬ 
bination of LK cycles, tidal friction and the expansion of 
the stellar radius caused by stellar evolution. This calcu¬ 
lation suggests that the number of BSSs formed through 
a combination of LK dynamics and stellar evolution in 
our models is negligible. 

Finally, note that we did not consider primordial triples 
that could somewhat increase the number of mergers 
from stellar triples shown in Figure 1111 However, we 
think that mergers occurring in primordial triples will 
have a negligible contribution to BSSs compared to merg¬ 
ers from dynamically assembled triples. In fact, due to 
the high stellar densities of GGs the outer binaries in pri¬ 
mordial triples are likely to be involved in several dynam¬ 
ical encounters which will lead to their disruption in the 
first few Gyrs of c luster evolution. Th e majority of triples 
in the catalog of iRiddle et al.l (|20I5ll have outer period 
P 2 > 10“^ days, which according to Equation (TT^ corre¬ 
sponds to a typical survival time Tenc ^ 10® yr within 
the core of a GC a nd 5. 10^° in the clus ter outskirts 
(see also Fig. 2 of iPerets fc Fabrvckvll2009ll . Moreover, 
given the old age of GGs stellar merger products that oc¬ 
curred in primordial triples would have already evolved 
off the main sequence, so that only dynamically formed 
triples could then produce BSSs at later stages of cluster 
evolution. 

7. SUMMARY 

In this paper we studied the long-term evolution of 
dynamically assembled triples in GGs. The initial condi¬ 
tions for the triples were obtained directly from detailed 
Monte Carlo models of the clusters. These systems con¬ 
sisted mostly of triples containing an inner stellar binary 
and triples containing a pair of BHs. The triple initial 
conditions were integrated forward in time using a high 
accuracy three-body integrator which incorporates rela¬ 
tivistic corrections and tidal terms to the equations of 
motion. The direct integrations allowed us, for the first 
time, to put constraints on the contribution of the LK 
mechanism to the population of coalescing BH binaries 
and BSSs in GGs. The main results of our study are 
briefly summarized below. 

(1) We find that for the majority of dynamically assem¬ 
bled triple systems in GGs, the time-scale for changes of 
angular momentum during the high eccentricity phase of 
a LK cycle becomes shorter than the orbital time-scales. 
We conclude that the standard orbit average equations 
of motion, often employed in previous studies, cannot 
accurately describe the evolution of such triples. 

(2) BH triples assembled dynamically in GGs often 
evolve such that the inner binary angular momentum 
changes on timescales of order or shorter than the inner 
binary orbital period (see Figure [2]). This can lead to 
the formation of eccentric GW sources in the frequency 
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band of aLIGO detectors (see Figure ED- 

(3) We estimated a realistic aLIGO detection rate 
of BH binary mergers due to the LK mechanism of 
^ 1 yr“^. We also compared this rate to that for BH 
mergers due to other processes (see Table 1). We hnd 
that mergers induced by LK oscillations are likely to be 
the main contributor to in-cluster GW sources for aLIGO 
detectors, and represent ~ 1% of the total detection rate 
of BH binary mergers assembled in GGs. 

(4) About 20% of all LK mechanism induced BH merg¬ 
ers in GGs will have an eccentricity larger than ~ 0.1 
at 10 Hz frequency, such that eccentric waveform tem¬ 
plates and/or dedicated search strategies will be needed 
for their efficient detection (see Figure |4|). Most of these 
eccentric sources will have an extremely large eccentricity 
(1 — e < 10“^) at the moment they first enter the aLIGO 
frequency band. Assuming perfect waveform templates, 
we Hnd a possible detection rate for eccentric mergers as 
large as a few events a year with aLIGO. We also show 
that the triple channel is likely the dominant formation 
mechanism for the production of eccentric inspirals that 
will be potentially detectable by aLIGO. 

(5) In the case of triples containing a stellar binary, 
the inner binary angular momentum often changes on 
timescales shorter than the tertiary orbital period, but 
rarely shorter than the inner binary orbital period. This 
implies that “clean” collisions between two stars, in 
which all passages prior to the collision are greater than 


the tidal dissipation scale, are rare (see Figured]). 

(6) We estimated the number of BSSs expected to form 
through mergers induced by the LK mechanism as well 
as those formed through other processes, including bi¬ 
nary stellar evolution, binary-single scattering and di¬ 
rect collisions (see Figure dH)- We find that the con¬ 
tribution of LK induced mergers to the population of 
BSS populations in GGs is typically < 10% of the to¬ 
tal. Even in clusters with relatively low central densities 
(- I02-3 Mqpc- 3) only up to - 10% of BSS binaries 
could have formed in dynamically assembled triples. 
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